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SEASONAL VARIATIONS IN THE F 2 REGION 


H. G. Mayr 
K. Mehajan 


ABSTRACT 

Semiannual and annual variations in the F 2 -region are discussed based on 
foF 2 and radar back scatter observations during the descending phase of solar 
activity (1958-1965). The semiannual effect appears as a persistent feature of 
the ionosphere which, demonstrably, is not related to fluctuations in the 10 k 7 cm 
noise or the EUV radiation. This gives support to theories that attribute the 
semiannual effec to variations in the lower atmosphere. Although theoretically 
predicted temperature variations (Volland, 1969) could quantitatively account for 
the observed semiannual variations in the height of the F 2 -peak, the variations 
in N m F 2 require additionally significant variations in the neutral composition at 
lower heights. Evidence for this is found in rocketborne [0]/[0 2 ] measure- 
ments at 120 km which show maxima during equinox and a maximum to minimum 
ratio of two consistent with the ionospheric behavior. The annual variation in 
the ionosphere, showing a winter to summer enhancement in the F 2 peak density 
at low latitudes and during low solar activity, can be explained by meridional 
winds at ionospheric heights. These winds also have a significant effect on the 
height of the F 2 peak and appear to be effective throughout the solar c:y cle as 


iii 



evident from radar backscatter observations, which show the peak at higher 
altitudes in summer than in winter. This is consistent with the observations oi 
King et al., 1968, pointing out that the electron density in the topside ionosphere 
is persistently higher in the summer hemisphere than in the winter hemisphere 
thus suggesting that in the upper ionosphere the wind effect on the height of the 
F 2 maximum has masked the processes responsible for the winter anomaly 
in N. 
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SEASONAL VARIATIONS IN THE F 2 REGION 

INTRODUCTION 

There are two major anomalies in the seasonal behaviour of the daytime F- 
region. The first and the more widely discussed one is the winter anomaly where 
the winter F-region critical frequencies (foF2) are found to exceed those in sum- 
mer, in particular at the higher latitudes and during the high solar activity (see 
e.g. Ratcliffe and Weekes, 1960 and references therein). The second, but not as 
widely discussed, is the semi-annual effect in foF2, with maxima in March and 
October and minima in July and January (see e. g. Burkard, 1951, Bhargava, 

1959, Yonezawa and Arima 1959, Yonezawa, 1959). The winter anomaly has, in 
fact, been demonstrated to be really an electron density depression in summer, 
(see e, g. Wright, 1963) and has been proposed by Johnson (1964) and King (1964) 
to be due to molecular oxygen abundance in the summer hemisphere, arising from 
the global upper atmosphere circulation with air rising over the summer and 
settling over the winter hemisphere (Kellogg, 1961). While this argument has 
been re-supported by Duncan (1969) for the winter anomaly in his recent work, 
no attempts seem to have been made in the past to explain the semiannual effect. 
In this paper we shall re-study the seasonal variations in the daytime F-region 
and try to explain the semi-annual effect in foF2. We shall discuss the overall 
behaviour of the F-region in view of three kinds of mechanisms, namely: (i) 
neutral-ion drag effects due to horizontal meridional winds (King and Kohl, 1965), 
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(ii) changes in the neutral temperature (Volland 1969) and (iii) changes in the 
neutral composition, for which we shall present some experimental evidence. 


DISCUSSION 

The data employed in our investigation were bottomside sounder measure- 
ments of the critical frequency (foF 2 ) and radar backscatter observations of the 
height of the F 2 maximum (h^) at Arecibo. In particular, the monthly means of 
these two parameters are considered, thus eliminating short terra variations. 

We chose the data from 1600 local time which appear to be closest to the day- 
time maximum during the yearly variations. This choice was made to justify 
our steady state approach in solving the ion continuity equation for the F 2 region. 
However, this point is not critical because we have not attempted to reproduce 
actual values for height and density of the F 2 maximum but rather to match the 
yearly variations of these two parameters; and on this long time scale the as- 
sumption of steady state is certainly appropriate for the ionosphere. 

Semiannual effects are observed both in the neutral and ionized components 
of the atmosphere. Burkard (1951) discovered first the semiannual effect in the 
F 2 -region. Patzold and Zschorner (1960) detected it in the thermospheric 
density, and subsequently semiannual variations were observed in the geomag- 
netic S q variation (Wagner, 1968) in the geomagnetic activity (Priester and 
Cattani, 1962) in the height of the F 2 maximum (Becker, 1966), in the electric 
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and neutral components of the mesosphere (Lauter et al., 3 966) and in the wind 
circulation within the mesosphere (Kochanski, 1963) and the stratosphere (Quiroz 
and Miller, 1967). 

Several solar-terrestrial interaction mechanisms have been invoked to 
explain these effects: PUtzold and Tsehdrner (1960) and Priester and Cattani 
(1962) attribute the semiannual variation to the solar wind impinging on the 
magnetosphere. Chandra and Krishnamurthy (1968) pointed out that a correlation 
exists between the 10.7 cm solar radio noise and thus EUV radiation, and the 
semiannual variation in the gas temperature. Others invoke processes that oc- 
cur in the lower atmosphere: Newell (1968) attributes it to a change in the meso- 
spheric circulation system and Volland (1969) shows that the semiannual vari- 
ations in the thermospheric temperature can be explained as resulting from tidal 
waves that are triggered by the .esospheric wind system. It will be shown here 
that ionospheric observations are consistent with the latter group of theoretical 
models. 

Figure 1 reveals that a semiannual effect in N m was present throughout the 
entire descending phase of solar activity. It is most pronounced at low latitudes 
with a ratio between maximum to minimum density of the order of 1.6. The vari- 
ation is barely recognizable at midlatitudes. 

As an indicator of the solar extreme ultraviolet, the 10.7 cm flux is shown 
for comparison. Its fluctuations reveal a maximum to minimum ratio of less 
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than 1.2. To the extent that the 10.7 cm flux represents the variations in the 
EUV radiation, its variations should correspond to fluctuations in both the 
thermospheric gas temperature and the electron production rate. An increase 
in the temperature enhances N 2 and 0 2 more than O, thus tends to reduce the 
F 2 maximum density. However, the increase in the photoionizing EUV flux in- 
creases the electron density proportionally. Thus the electron density variations 
induced by fluctuations in the EUV (10.7 cm flux) are opposed by the increase in 
T .vnich accompanies the flux-increase, and one therefore expects a smaller 
variation of N m F 2 than 1.2. This value is significantly smaller than the observed 
ratio in the semiannual effect and we consider this as one evidence that the semi- 
annual effect cannot be attributed primarily to the solar EUV radiation. It is, 
however, clear that EUV fluctuations are not negligible when superimposed on 
the semiannual variations and this may be responsible for the positive correla- 
tion that sometimes exists between the 10.7 cm flux and the electron density in 
periods wher this flux shows a quasi-half year period. Additional evidence that 
the fluctuations in the 10.7 cm flux are not responsible for the semiannual effect 
(as a prevailing phenomenon) can be derived from the fact that summing the mean 
monthly 10.7 cm values over a solar cycle (1958-1969) does not reveal a signifi- 
cant semiannual effect. A similar conclusion was arrived at by Newell (1968). 

Volland’s (1969) theory of the semiannual effect in the thermosphere is 
ba^d on a wave energy input from the lower atmosphere that produces an almost 
constant temperature increase at all altitudes in the thermosphere. Figure 2 
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shows N m and h m calculated (see Appendix) for a 100° temperature enhancement 
for minimum and maximum solar activity conditions. Furthermore, we show 
these parameters as functions of the composition ratio [0)/[0,l at 120 km. 
Assuming this to reflect only variations in O we adopted a similar variation in 
the |0]/[N 2 j ratio. 

From Figure 2 it is apparent that a temperature increase of 100°K can only 
produce an electron density maximum to minimum ratio of 1.15, small compared 
with the observed value of 1.6. The temperature effect on h m is significant how- 
ever as it increases h m by 35 km during minimum solar activity and by 45 km 
during maximum solar activity. This appears to be consistent with observations 
by Becker (1966) which show semiannual variations in h m of this order. 

As the [ O ] / [0 2 ] and [O ] / [N 2 ] ratios significantly affect the F-region 
electron densities, we have collected all the neutral composition measurements 
to date (see Table 1) to look for any semiannual effect in these quantities. Figure 
3 shows a definite indication of a semiannual effect in the [O] / [0 2 ] ratio at 
120 km, although some of this variability may be attributed to the somewhat 
different measurement techniques employed in these rocket experiments. 
Furthermore, neutral composition measurements (References from Table 1) 
show significant variability in the [ O ] / [N 2 ] ratio (not shown here) while 
simultaneously [ N 2 ] remains constant. Unfortunately, the few data available 
on the latter ratio are insufficient to provide information on its semiannual 
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component* Variations in the neutral composition as described in Figure 3 are 
more than sufficient to produce the semiannual variations in the F 2 -maximum 
density as demonstrated in Figure 2). The observed variations in the electron 
density are in fact less than the composition amplitude would require. 

Annual Variations 

During low solar activity and at low latitudes N m is larger in summer than 
in winter (Figure 1). Furthermore, radio backscatter observations of h m in 
Puerto Rico (Figure 4) show that this parameter increases by about 100 km from 
winter to summer. 

Composition changes as shown in Figure 3 cannot account for these effects 
because they are not in phase with the summer enhancement and, moreover, they 
do not affect h m appreciably. Furthermore, reasonable variations in the gas 
temperature would also not produce enough variation of both electron density 
andh„,. 

King and Kohl (1965) and many other investigators have stressed that hori- 
rontal meridional winds in the neutral atmosphere significantly affect the distri- 
bution of the ionospheric plasma. Kohl and King (1967), Geisler (1966) and 
Volland and Mayr (1969) deduced the wind field from the empirical atmospheric 
model that Jacchia (1964) derived from satellite drag measurements. A gross 
feature of this wind field is that the winds blow away from the atmospheric pres- 
sure bulge that forms at the subsolar point (Jacchia and Slowey, 1966) and that 
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the wind velocities increase away from the subsolar point toward the poles. 
During the year the subsolar latitude moves between 23.5 north and south geo- 
giaphic. At American longitudes (where the ionosphere data were taken) the 
subsolar point is at 37 C N geomagnetic latitude at summer solstice and at 11 # S 

at winter solstice. Thus at low latitudes in summer the daytime winds blow 
toward the equator. As only the wind component parallel to the magnetic field 
affects the ionosphere, these equatorward winds produce at low latitudes an up- 
ward drift of ionization that enhances height and density of the F 2 peak. Con- 
versely, during winter the winds blow away from the equator and produce corre- 
sponding downward ionization drifts that decrease both height and density of the 
F 2 -maximum. Thus winds can reasonably be invoked to produce winter to 
summe*' enhancements in the ionization and peak height of the F 2 -region. 

The approach of Volland and Mayr (1969) is here adopted to describe the 
structure of the meridional thermosphere winds by means of spherical functions 
of low degree: 


21 


u = U 10 P io cos X + U 2i v cos (co(T - r )) 

COS A. 


(1) 


where 

P 10 = sin \ 

p 21 =| sin (2A) 
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K - geographic latitude 


s 

t - r + — (local time) 

(A) 

t = universal time 

O) (angular frequency with r Q the period of one day) 

T o 

i - geographic longitude 

r 21 = 1400 local time (observed diurnal maximum in atmospheric density) 

u is positive in fourth direction. 

Equation (1) is only valid for equinox conditions but with an appropriate trans- 
formation it can be formally adapted to describe also winter or summer condi- 
tions. In this paper we shall not go into the details of this modification, but 
rather restrict ourselves to the very simple case where the diurnal density 
maximum occurs at noon geomagnetic local time (at American longitudes) . Then 
we are justified to replace X. by (k - 77 ), where 77 is the inclination of the sun to 
the geographic equator (77 = 23.5° for summer and 77 = -23.5° for winter). At 
American longitudes the magnetic equator is south of the geographic equator and 
thus we can describe the wind field by substituting in equation (1) k = {$ + 13°) 
where we define 6 on the geomagnetic latitude. With I the magnetic dip angle, 
the wind component parallel to the magnetic field is u ( | = u cos I. The vertical 
wind component, that enters into the ion continuity equation is then 

u l = U H s ^ n * = u cos * s ^ n I » 
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Its latitudinal distribution is shown in Figure 5, for summer, winter find equinox 
conditions in the northern hemisphere. 

For illustrative purpose, we show in Figure 6 several electron density 
profiles which were computed for minimum solar activity conditions (at 30° geo- 
magnetic latitw e) employing upward and downward wind components. From this 
diagram it if-, evident that the winds, shown in Figure 5, can quantitatively account 
for the observed variations in h max (Figure 4) and N max (Figure 1) during low 
solar activity. A feature of the electron density distribution in Figure 6 is that 
upward winds are much more effective in increasing N m and h m than downward 
winds are in decreasing these parameters. (This is to be expected since 
downward winds bring the F 2 -layer closer to the chemical equilibrium region.) 

For this reason it is understandable that the wind effects are most pronounced 
in N m at low latitudes where, according to Figure 5 the upward wind components 
are largest. 

Our interpretation of the summer enhancement of N m by atmospheric winds 
appears only satisfactory for minimum solar activity conditions. Figure 1 reveals 
that during high solar activity the semiannual effect prevails at low latitudes 
while the winter anomaly, with higher density values in winter than in summer, 
prevails at midlatitudes. 

It was pointed out by Geisler (1966) and Volland (1969) that the magnitude 
of the winds are controlled primarily by ion-neutral drag rather than viscosity 
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affects. Thus, at high solar activity the greater electron density tends to in- 
crease the ion drag and to decrease the wind velocity. But at high solar activity 
the pressure in the neutral atmosphere is also enhanced and this would tend to 
increase the wind effect in F 2 region heights. The observations of the height of 
maximum (Figure 4) show that the ratio of the winter to summer enhancement 
did not change significantly between 1958 and 1965, which suggests that wind 
effects are of comparable importance throughout the solar cycle. Therefore, 
other mechanisms must be invoked at high solar activity to selectively depress 
and mask the summer enhancement of N m (but not h m ) that is so apparent during 
low solar activity. 

As shown earlier, variations in the neutral composition can affect N m without 
considerably changing h m . Thus, we infer that during high solar activity [0]/[0 2 ] 
and [0]/[N 2 ] decrease from winter to summer* and thus decrease N m in summer 
relative to its value in winter. This effect could counter balance or even neu- 
tralize the summer enhancement that would normally be produced by wind effects. 
Therefore it is understandable that at low latitudes the semiannual effect pre- 
vails. At midlatitudes the wind effects are smaller and therefore the composition 
changes prevail with the result that the winter values of N m are enhanced over 
those in summer. 


Most of the measurements in ‘’'able 1 and Fig. 3 correspond to periods of low solar activity and 
thus we cannot check on this postulation. 


10 



The atmospheric composition changes required to produce the winter 
anomaly, are predicted in the circulation theory by Kellogg (1961), Johnson 

(1964a) and King (1964). These authors propose that solar heating in the lower 
ionosphere induces upwelling of atmospheric molecular species over the summer, 
and settling of atomic oxygen over the winter hemisphere. As the ionospheric 
heating increases with increasing solar activity, it is understandable that this 
mechanism would be primarily effective during maximum solar activity which is 
consistent with the existence of the winter anomaly. 

King et al. (1968) reported that topside sounder observations do not show the 
winter anomaly. This behavior supports further the wind effects discussed be- 
fore if we consider that the electron density in the topside ionosphere not only 
depends on N m but very significantly also on h m (see Figure 6). The latter of both 
reveals the wind effect throughout the entire period of changing solar activity 
(see Figure 4), thus tends to produce winter to summer enhancements in the 
electron density of the topside ionosphere. 

SUMMARY AND CONCLUSION 

We have discussed foF 2 observations over a wide range of latitudes and 
throughout the descending phase of solar activity (1958-1965), and radio back- 
scatter observations of the height of the F 2 maximum and Arecibo between 1958 
and 1965. 
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The semiannual variation in N m appears as a persistent feature of the iono- 
sphere. Fluctuations of solar EUV indicated by the 10,7 cm radio noise are 
found to be quantitatively insufficient to account for the observed amplitude of 
the semiannual N m variation, furthermore we do not see a prevailing semiannual 
component in the 10.7 cm flux over a solar cycle. For these reasons solar EUV 
must be excluded as the primary source for the semiannual effect. 

The ionospheric evidence supports theories that attribute the semiannual 
effect to variations in the lower atmosphere (Newell, 1968, Volland, 1969), To 
explain the ionospheric density variations, semiannual variations in the neutral 
composition are required. We find some evidence from the few available neutral 
composition measurements that the ratio [0]/[0 2 ] varies by a factor of two and 
is greatest during equinox. The phase and amplitude is consistent with the semi- 
annual variation in N m . 

m 

In a subsequent paper, an attempt is made to explain this composition effect 
by invoking the global wind circulation within the lower atmosphere (Mayr, 
Newton, Volland, and Mahajan, 1969). 

The annual component in the ionosphere showing a winter to summer en- 
hancement in N m at low latitudes and low solar activity, can be explained as re- 
sulting from horizontal atmospheric winds. Observations of the height of the F 2 
maximum at Arecibo show a strong annual variation at all levels of solar activity. 
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These annual changes in h m constitute important additional evidence for atmo- 
spheric winds, thus indicating that winds are effective throughout the solar cycle. 

The latter we suggest as an explanation for the fact that the winter anomaly 
is not observed in the topside ionosphere (King, et al., 1968). 
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APPENDIX 


It. is here assumed that 0 + is the major ion constituent within a neutral 
atmosphere where atomic oxygen, 0, is the predominant neutral constituent. The 
oxygen ion density is then conti’olled by photoionization of 0 with a production 
rate Q, charge transfer to molecule N 2 and 0 2 with a loss rate L[0 + ] , and dif- 
fusion of 0 + through 0 with a transport velocity V Q parallel to the magnetic 
field. In steady state, the continuity equation was the form 


0- L[0*1 - B-2- 

os 



0 ) 


where 

s is the distance along a field line and B is proportional to the intensity of the 
magnetic field, 

B = J_ (1 +3 cos 2 6) l/2 
r 2 

therefore accounting for the divergence of the magnetic field tubes. 

Conservation of momentum requires that the momentum equation 

e o*o to] t0 + ] (V Qt - V„) = - kT li°ll - mg,, ion (3) 

f* s 

is satisfied, where 

t?o+o = 2 x l<r 32 (cgs), drag coefficient for collisions between 0 + and o» 
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~ velocity component of 0 parallel to the magnetic field, 

T - (T e + T ),2 

T v ~ electron temperature, T t - ion temperature 
k = Boltzmann constant 
m - mass of 0 + 

g I( “ gravitational acceleration parallel to the magnetic field. 

Jacchia’s model (1965) was adopted for atomic oxygen, molecular nitrogen 
and oxygen, and gas temperature. The loss coefficient for charge transfer from 
y to . v 2 was chosen to be 2 x 10“ 12 cm 3 sec"* 1 and for 0 + to 0 2 to be 1 x 10" 11 
cm 3 see -1 . The photoionization probability was assumed to be 5 x 10“ 7 x F/80, 
where F is the 10.7 cm flux (10' 22 w/m 2 c/s). Attenuation of the ionizing EUV 
radiation has been rigorously taken into account. For the electron temperature 
we used 1500°K which is consistent with satellite measurements by Brace et al. 
(1965). The ion temperature was assumed equal to the gas temperatuie an as- 
sumption that is justified from energetic considerations. 

Integrating Equation (1) along a field line up to the equator yields 

V l0+ ] = §■ V, [0 + ] - B f ' (0 - L [0*] ) (4) 

e ' 

where subscript e refers to the field line intersection in the equatorial plane. 
Combining Equations (3) and (4) one arrives at the differential-integral equation 
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0*)V ll + B j" ' «?-LIO + ])^ - 

•'s 

Representing the right hand side of Equation (5) by D(Sj I O + ] ) and assuming this 
function to be known, this equation can then be integrated thus yielding the integral 
equation 


| Veto'll- (5) 

e 


kT 


ftf] 

ds 


+ mg,, [0 + ] 


o + o ‘ 0 ' 




icn 


[0 + ] 0 exp 




"D( Si [0 + ]) 

' [0 + ] ' 



(0) 


Equation (G) is solved for the following boundary conditions: 

1. At low altitudes 0 + is assumed to approach the value determined by 
photochemical equilibrium 



2, In case of symmetry with respect to the equatorial plane (equinox condi- 
tion) it is required that 



0 . 


The solution of Equation (G) was accomplished by means of an iterative 
procedure in which the function D(s x [0 + ])/[0 + ] was treated as a perturbation 
in the following scheme: 


For the k th solution the perturbation function was assumed to be 
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D(s, (0*l k ) DfSjlO*],.,) 

l0 \ H ' O ’ 1 .- 1 

where • Is a parameter that varied between 

0 s Vi ' e k < • £ 1 

For the initial solution we assumed * x - 0 which is equivalent with diffusive 
equilibrium. In subsequent solutions k is gradually increased to allow the per- 
turbation function to become effective slowly. This procedure assures the ap- 
plicability of the perturbation method which demands small effects from the 
"perturbing” function. A successive application of this procedure leads then to 
a convergence on the solution of Equation (6) when e k approaches 1. 

One can immediately verify that this solution procedure requires the more 
iterations the larger the chemical terms are when compared to the transport 
terms in the continuity equation. This imposes a lower altitude limit for 
reasonably fast convergence. Employing an. IBM 360-75 computer we started 
at 180 km and .ve achieved computer solutions accurate to better than 1% in 20 
seconds of computing time. 
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Table 1 

O and 0 2 Number Density Ratios at 120 km Measured by 
Rocket- Borne Mass Spectrometers 


V? 

Date 

Time 

Latitude 

n(0)/n(0 2 ) 

Reference 

175 

Sept. 23, 1960 

0050 LT 

'Middle' 

0.72 

Pokhunkov (1963a, b, c) 

95 

May 18, 1962 

1302 EST 

38° N 

(1.2) 

Schaefer and Nicholas (1964) 

73 

Mar. 28, 1963 

0235 LT 

38° N 

1.2 

Schaefer (1966) 

77 

June 6, 1963 

0730 MST 

33 °N 

0.78 

Nier et. al. (1964) 

82 

Nov. 26, 1963 

1316 LT 

38°N 

1.2 

Schaefer (1967) 

72 

Feb. 18, 1965 

1409 LT 

59° N 

0.75 

Schaefer (1967) 

72 

Feb. 19, 1965 

0317 LT 

59° N 

0.75 

Schaefer (1967) 

75 

Apr. 15, 1965 

0345 MST 

33 °N 

0.33 

Hedin and Nier (1966) 

75 

Nov. 30, 1966 

0445 MST 

33°N 

(0.5) 

Kasprzak et. al. (1968) 

98 

Dec. 2, 1966 

1409 MST 

33° N 

0.46 

Kasprzak et. al. (1968) 

76 

Dec. 11, 1965 

0505 MET 

40°N 

1.56 

Mauersberger et. al. (1967) 

163 

Dec. 12, 1966 

1320 CST 

59° N 

0.87 

Gross et. al. (1968) 

119 

June 21, 1967 

1249 MST 

33°N 

0.55 

Krankowsky et. al. (1968) 

131 

July 20, 1967 

0200 MST 

33°N 

0.41 

Krankowsky et. al. (1968) 

131 

July 20, 1967 

1224 EST 

33°N 

0.41 

Krankowsky et. al. (1 968) 
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Figure 2. Computed N and h as functions of the l0}/l0 o ] ratio at 120 km varying from 
one to two, and furthermore depending on a constant temperature increase in the entire thermo- 
sphere (T is the exospheric gas temperature). Solid lines show the variability during minimum 
and maximum solar activity. The thin lines demonstrate the variations during the increasing 
phase of solar activity. 




MASS SPECTROMETRIC MEASUREMENTS AT 120 km 

(See Table 1) 
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Figure 3- [0 ]/ [Oj] ratio at 120 km plotted from measurements referred to in Table 1. Note that the majority of 

data points exhibits a few annual variation (observed live). 







SUMMER 



(oas/iu) T n 


29 


Figure 5. Vertical component of the meridional wind field shown for 
equinox, summer and winter in the northern hemisphere. 









